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ABSTRACT 


The existence of a solution defined for all t and possessing 

» 

a type of boundedness property is established for the perturbed non- 
linear system y = f(t, y) + F(t, y) . The unperturbed system 
x = f(t, x) has a dichotomy in which some solutions exists and are 
well behaved as t increases to « and some solution exists and 
are well behaved as t decreases to A similar study is made 

for a perturbed nonlinear differential equation defined on a half 
line, say R + , and the existence of a family of solutions with 
special boundedness properties is established. Finally, the. ideas 
are- applied to the study of integral manifolds. Examples are 
given. 
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1. INTRODUCTION 

'The following is a study of the system 

(1) y = f(t,y) + F(t,y), 

which is regarded as a perturbation of the nonlinear system 

(2) x=f(t,x). 

We impose hypotheses on (2) which guarantee the existence of a 
bounded solution (or a family of bounded solutions) and prove that, 
under conditions on the perturbation terms, such solutions are also 
present in (l). Results of this type have been obtained by May [11] 
and Harbertson and Struble [8] for nonlinear systems and Coppel [3], 
Hallum [ 6] , [7] and Hale [ 4] , [5] for perturbations of linear 
systems. The present work removes certain hypotheses from the past 
results and introduces function spaces which allow new types of be- 
havior to be studied. 

Theorem 1 concerns the case when (l) and (2) are defined 
for all t and give the existence of particular solutions defined 
on R. Theorem 2 covers the case when the systems are defined on 
some half line t > a and concern the presence of special solutions 
on this interval. Theorem 3 extends the above ideas to integral 
manifolds. Examples are discussed. 
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• 2. NOTATION AND GENERAL SETTING 

Let p and q be nonnegative integers with p + q = m > 0, 
let I be an interval of the type I = {t > a} (a may be -»), 
let D^(t), Dg(t) be continuous nonsingular p X p and q x q 
matrices respectively on I and for 0 > 0 define 


ft^ = {(t, x) in I X R P : |D 1 (t)x| < 0 ), 


We assume 


l ftp = {(t, x) in I X R q : |D ? (t)x| < 0 } . 


that (2) may be written in the form 


( 3 ) 


*1 = 

x 2 = f 2 (^ x 2)> 


where for some 0 > 0, f^ is a continuous function on ft^ into 

R p , fp is a continuous function on ftp into R^, where f^ and 

fp have continuous partial derivatives in x^ and Xp 

respectively and f^(t,0) = 0, fp(t,0) = 0. For (x, a^) in ft^, 

(x, ap) in ftp we denote the solutions of (3) such that 

x^(t) = a^, Xp ( t) = ap by x^(t, x,a^), Xp(t, x,^). We assume that 

for some 0 < r < o and (x^a^) in ft' in ftp, x^t, x, a^) 

exists for t > t and Xp(t, x,ap) exists for a < t < t. 

£ 

For convenience when a is in R 


we will use norm 
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*th 

| a| = max i\\\) where i = 1 , 2, . . ., 4, i component of 

i 

a. Then for in R , a 2 = R^, a = col(apa 2 ) 

| a| = max { | aj , | ag| } . Let he the continuous functions y on 

I into R^ with |D^(t)y(t)| < y and B 2 he the continuous func- 
tions y on I into R^ with | D^C^h) y('t) | < Y. 

bx. 

We denote the matrix (t, T, a.) hy ®^(t, T, a.) , i =-1,2, 

i 

and assume that F^ and are continuous functions on 

n r = {(t,x 1 ,x 2 ) in I X R m : |D 1 (t)x 1 | < r, |D 2 (t)x 2 | < r} 

into R 15 and R q respectively. 

3. SYSTEMS DEFINED OR R 
section we consider the case when I = R. We note 

Y* y 

system (3) has a solution^in B^ X B 2 , viz. 

We seek hypotheses on the perturbed system 

y-L = ^(^y^ + F 1 (t.y p y 2 ), 

(4) 

y 2 = f 2 (t,y 2 ) + F 2 (t,y p y 2 ), 

Y y 

which guarantee the existence of a solution (y^y^) in B| x B^ . 


In this 
the unperturbed 

X 1 = °> x 2 = °* 


Theorem 1„ Assume that we have 
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t 

/ •|P 1 (t)® 1 (t, s,y 1 (s))F 1 (s,y 1 (s),y 2 (s))| ds < r, 

-00 . ' 

.(5) * . 

/ |D 2 (t)® 2 (t,s,y 2 (s))F 2 (s,y 1 (s),y 2 (s))|ds < r. 


V Y 

for all t in R and (y^y^) in X B^. There is a solution 
(y^^y^t)) of (4) defined for. all t and (y x ,y 2 ) is in 


Proof : For each positive integer k, let he the set of func- 

tions (y^ y 2 ) where y^ and y p map [ - k, k] continuously into 
IT and respectively. For y = (y^ y^) in let 


I y| = max {sup |D (t)y (t)|, sup | D ? (t)y p (t)| }, 
t t 

then (B^, ||) is a Banach space and B^ X B p is a closed convex 
subset. (The functions in B| X B^ are here restricted to [-k, k].) 

y y 

On B^ x B p we define the transformation T by Ty = u where 


( 6 ) 


u^t) = / ® 1 (t,s,y (s))F 1 (s,y (s),y 2 (s))ds, 
-k 


^(t) = -/ ® 2 (t,s,y 2 (s))F 2 (s,y 1 (s),y 2 (s))ds, 
t 


for -k < t < k. It is an easy exercise to check that Schauder's 

r y 

fixed point applies so there is a function y in B^ X B p 
defined for -k < t < k with Ty = y. We have 
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y L (t) = F 1 (t,y 1 (t),y 2 (t)) + / s)F 1 (s,y 1 (s),y 2 (s))ds, 


-k 

t 


y 2 (t) = F 2 (t,y 1 (t),y 2 (t)) + / H 2 (t,s)F 2 (s,y 1 (s),y 2 (s))ds, 

k 


where 


df. 


H^s) = - (t,x i (t,s,y i (s))® i (t,s,y i (s)). 


i = 1,2. Here we have used the well known ([ 3 ], page 22) theory of 
the varialional equation. Also 


f 1 (t,y 1 (t)) =/ [f 1 (t,x 1 (t, s,y 1 (s)))]ds, 

“k 


= f H x (t, s)[y 1 Cs) - f 1 (s,y 1 (s))]dsj 
-k 


and a similar expression holds for f 2 (t, y 2 ^)) • For i = 1*2 and 


v.(t) = y ± (t) - f^y^t)) - F i (t,y 1 (t),y 2 (t)), 


we obtain 


w,(t) = / H..(t, s)w n (s)ds, 
-k 


w 9 (t) = / H 9 (t, s)w (s)ds, 


(7) 
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■which implies 'w^(t) = 0, w^t) = 0. Consequently y = (y^y^) is 
a solution of (4), -R < t < k in X Bg. 

00 

Let [y(t, k)} be a sequence of fixed points of T on 
X— X 

Y* Y* 

x * ^2 ( res- tricted to [-k, k]). The following statements show 
there is a function y in B^ X B^ and a subsequence [y(t, n^) }^_-l 
such that 

(8) lim y(t,n^) = y(t), 

, k -» <» 

and the limit is uniform on compact t intervals. 

There is a subsequence [y(t,n.^)} converging uniformly on 
[-1,1] since the original sequence is uniformly bounded and each 

function satisfies (4) on this interval. Similarly, there is a 

•> 

subsequence [y(t, n 2k ) } “ =1 of the sequence [y(t, n ) }” =1 con- 
verging uniformly on [-2,2], In this way we obtain a chain of 
subsequences (y(t, n converging uniformly on [-j,j]. Put 
y(t, n^) = y("k>n^) to obtain (8). Since each y(t,n^) is a solution 
of (4) so is y. 

May [11] gives a similar theorem when D^(t) = I, DgCt) = I; 
however he requires the additional hypothesis of a Lipschitz condition 
in z on $(t, s, z)F(s, z) . The essential difference in the proofs is 
that we define T using a finite limit k and use functions on 
intervals [-k, k]. The proof that the fixed point y satisfies (4) 
reduces to the fact (7) implies w^w^ = 0, an easy observation when 


k is finite. 
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Hallum [ 6], [7] and Hartman and Onuchic [9] introduce the . 
matrix D(t) = g(t)l in studies of a perturbed linear system on a 
half line. Here g(t) is a continuous nonnegative functions. 

Condition (5) is restrictive on the linear parts in x of 
the fp fg functions. Suppose, for example, f^tjX^) = A^C^x^ 

+ h(t,x^) where h(t,x^) = o(| xj ) uniformly in t as x^ -» 0. 

Let X^(t) be a fundamental matrix of solutions for x^ = A^(t)x^, 
then since <J>(t, s,0) is a solution of z = f (t,0)z, z(t) = I, 
we see (4) implies, in particular, 

t . , 

/ |D 1 (t)X 1 (t)X 1 x (s)F 1 (s,0,0)|ds < r. 

-00 

t _1 

(Theory concerning systems with hypotheses f |X(t)X (s)|ds < K 

-00 ' 

is given by Coppel [3].) If, in addition, A^(t) =0 we see 
F^(s,0,0) must be integrable. 

As an example take p = 1, q = 0 and consider the system 

c, 

00 

(9) y = a(t)y - Z b.(t)y 2l+1 + F(t,y), 

i=l 

where a(t) is a continuous function on R to be further restricted 

later, the b^(t), i = 1 , 2 ,..., are continuous nonnegative functions 

on R, the series Z b^(t)y converges to a continuous function 

2 ± 

with continuous partial derivative in y given by Z b.(t)y , 

2 

and the function F is continuous on R and will be further 
specified later. We take as the unperturbed system, 
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W 4 

x(t) = a(t)x - E b.(t)x 1+ , 
i=l 1 


then the variation equation is 


z = [a(t) - E h i (t)x 2i (t) ]z. 


Consequently, 

t 

/ a(s)ds 

‘ |®(t,T,r)| < e T 

Case 1 . Suppose a(t) = -1 (or any negative constant). 

Then for D(t) = 1 we have that if there is a a > 0 such that for 
continuous functions y(t) with |y(t)| < a, 

/ e"^ t " S ^|F(s,y(s))| ds < a, 

-00 

then (9) has a solution y* defined for all t with |y*(t)| < a. 

We notice here the nonlinearities E h^(t)y 2 ^ + ''' are "harmless". 

This observation shows the nonlinear theory allows much larger bounds 
for some systems than the corresponding linear perturbation theory 
gives (Coppel [3], page 137). It is interesting to note that for 
a(t) as given in Coppel [3], page 73, the above conclusions also 
hold. In this case the basic linear system is not exponentially 


stable. 
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Case 2 . 
D(t) .function. 


Sometimes the form of the linear term dictates the 
Suppose a(t) is given by 


a (t) 


f 


2t 


t 2 + 1 


< 


-1 

t + 2 


t < 'I, 


-1 < t. 


Then it is easy to see there is a k > 0 so that for 

I t 

D(t) = k/(,|t| + 2), J D(t)$>(t, s, y(s))ds is bounded, say by M, 

| -00 

for all t and continuous functions y. If there is a a > 0. 

So y in implies |F(t,y(t)| < a / M then ,(8) has a solution 

y* defined for all t with |y*(t) | ^ ct(| t| + 2) /k. Other 
choices for a (t) will give decreasing D(t) fractions as t -* °°. 


Case The form of the perturbation term may dictate a 

D(t) fraction. Suppose a(t) = -1 and F(t, y) = h(t)k(t, y) where 



t < 0, 
t > 0. 


Then 
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/ e ^ S ^h(s)ds = 




1 . -t 

2 + t e , 


t < 0, 
t > 0, 


consequently a natural choice for D(t) is 


t < 0, 
t > 0. 

In most instances the solutions of (3) and consequently 

I 

<I>p $2 are not known precisely. However, as in the example above, it 
may be possible to obtain information which implies the hypothesis 
of Theorem 1. Consider the following situation in which, for con- 
venience, we assume q = 0, and suppress the subscript 1 notation. 
Suppose D(t) is given and there are positive numbers a,K, a set 
!1CE^ and a continuous real valued function X on R such that: 

(a) for s in R, | t| < | D 1 (s) | a, x(t, s, r) lies 
in ft for t > s; 

(b) p[ f x (t,x(t)) ] < X(t) for all continuous functions 
■ x from R to ftj 

t t 

(c) D(t) / exp / \(u)du ds < K. 

-00 s 

We notice that for y in B°, | y( s) | < | D s) | a, hence 


2e 


D(t) = 


1 + 2te 


-t * 
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x(t, s,y(s)) is in H for t > s. By a well known result . (Coppel 
[3.]j page 58) condition (b) implies 

t 

|®(t,s,y(s))| < exp / n(f x [u,.x(u, s,y(s))])du; 

s 

hence condition (c) together with a boundedness assumption on F 
will imply inequality (5) in Theorem 1. 

4. SYSTEMS DEFINED ON A HALF LINE 

In this section we consider the case a finite. We note 

the unperturbed system may have a family of bounded solutions on 

t > a. We seek hypotheses on the perturbed system (4) which 

T Y 

guarantees the existence of a family of solutions in X B 2 . 

Theorem 2 . Assume that we have 

t 

/ |D 1 (t)<I> 1 (t, s,y 1 (s))F 1 (s,y 1 (s),y 2 (s))| ds < r/2, 

( 10 ) 

00 

/ |D 2 (t)d> 2 (t, s,y 2 (s))F 2 (s,y 1 (s),y 2 (s))| ds < r, 

*fc 

Y* v* 

for all t in R and (ypY 2 ) in B^ X Bg. Further assume that 
for some A > 0 and | a^| < A, a^ in R^ we have 

D 1 (t)x 1 (t,a,a 1 )| < r / 2 . 


( 11 ) 
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Then for | aj < A , in R^ there is a solution (y^(t),y 2 (t)) 

of (4) defined for ^t > a, (y-|_,y 2 ) is in B^ X B 2 and y^a) = a.^. 

Proof : For each positive integer k > a, let B^ be the set of 

functions (y^yg) where y^ and y 2 map [ a, k] continuously 
into R^ and R^ respectively. On B^ X B 2 restricted to [a, k] 
we define the transformation T by Ty = u where 


t 

^(t) = x 1 (t,a,a 1 ) + / 3> 1 (t, s,y 1 (s))F 1 (s,y 1 (s),y 2 (s))ds, 



t) 


k 

-/ 1> 2 (t, s,y 2 (s))F 2 (s,y 1 (s),y 2 (s))ds, 
t 


for a < t < k. As before such a transformation has a fixed point 

which is a solution of (4). The family of such' fixed points 

» 

(a < k < “) has a convergent subfamily which converges to a func- 
tion (ypy 2 ) satisfying the conclusions of the theorem. 

Suppose that, in addition to line 2 of inequality (10), we 
have a system such that 


t 

(12) / |B 1 (t)«D 1 (t, s,y 1 (s))| ds < Y* 

0 c 


for all y^ in B^\ Then appropriate boundedness conditions on F^ 
will give line 1 of inequality (10). Inequality (12) for D^(t) = I 
gives a type of exponential stability to x(t,Q£p a^) if 
ff(t, Xf) = A^(t)x^. Coppel [3], page 68, proves that there is a 
constant R such that 
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7?(t-a) 

tx 1 (t)x“ x (a)| < Ne , a+l<t. 

In this case inequality (ll) is not an added restriction hut merely 
defines notation. Corresponding theorems are given for general 
D(t) by modifying Hallum [6], page 255-6. Brauer [1], gives a 
theorem in the direction for nonlinear functions f^(t,x^) which 
could again be modified to include D(t). 


5. INTEGRAL MANIFOLDS 

, r 

In this section we consider the system 

e = h(6,t,z) + H(e,t,y,z, e), 
(13) y ='f(t,y) + F(0,t,y,a,e), 

z = eg(z) + eG(0,t,y,z,e), 


1 m n 

where (6, y, z) is in R X R X R and where hypotheses will be 
introduced to insure the existence, for small €, of an integral 
manifold of solutions. The form of the system and the hypotheses 
given' are motivated by previous work by Hale [4], [5] and Harbertson 
and Struble [ 8] . Such problems arise in the "method of averaging" 
introduced by Kryloff and Bogoliubov [10] (see also [2]) and studied 
extensively by many. 


The ideas introduced in the previous sections of this paper 


are applied to the study of ( 15) . A treatment of such a system 
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without the D matrix has been given in [8]; however our treatment 

improves the allowable bounds for the perturbation. In addition, 

> 

the form of the 0 equation has been changed to allow a larger class 
of examples. Such an example is given at the end of this section. 

It is possible to present this theory for the case where f 
is split into two functions, f = column ( f^, fg) as was done in the 
preceding sections. Correspondingly the g function (i.e., the z 
equation) can also be split into two pieces, one which is well 
behaved as t -» the other well behaved as t -» (see [8]). In 
order to present these ideas without unnecessary clutter we will not 
make these decompositions of the y and z equations. It will be 
clear how the hypotheses must be altered to obtain a corresponding 
theory with the dichotomies present in the y and z equations. 

Let D(t) and E(t) be continuous nonsingular mx m 
and n x n matrices respectively on R and for a = (a^, a^) in 
R + x R + we define 

n a = { (t,y, z) in R 1+m+n : |D(t)y| < |E(t)z| < a 2 ), 

= ((t,y) in R 1+m : | D(t)y) < 0 ^, 

= ((t,z) in R 1+n : J E ( t) z| < a 2 ) , 

and consider the following hypotheses. Let a in R x R and 
e Q > 0 be given. 
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(1) h is a continuous function from R^ X into R^. 

0 / (J 

H, F and G are continuous functions from R x ft x [0,6 ] into 
R^, R m and R n respectively, f is continuous on has a con- 
tinuous derivative in y and f(t, 0) =0. g is continuous on R n 

and has continuous derivatives in z. There is an cd = (cd, ,o> 01 . . . ,a> ) 
in R m such that h, H, F and G have period ox in 0^. 

H, F, and G vanish for 6=0 and |h| <C , |H| < Cg on their 

domains. ■ ' 

(2) The solutions of y = f(t,y), y(t) = a, (f,a) in 
are denoted "by y(t, t, a) and we assume for some 0 < a* < c , 

y(t, T,a) exists for t > t when (t, a) is in . We denote by 
$(t, t, a) the matrix ^ • 

(5) The solution of z = g(z), z(0) = b, b in R n . is 

denoted by z(t, b) and we assume z(t,b) exists for t > 0. We 

Sz 

denote by A(t,b) the matrix . 

(4) There is an N > 0 such that on its domain 

|h(e,t,z) - h(e*,t,z*)| < u(| z-z*| + | e-0*|). 

There is a continuous function jj(t, e) on R x [0,6 ] decreasing 
to 0 as e -» 0 such that on its domain 

|H(0,t,y,z, e) - H(8*,t,y*,z*, e)| 

< i(t, e) (| 0-0*| + |y-y*| + |z-z*|). 
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There is a continuous function V Q (e) on [0, e Q ], V q (o) = 0 and 
a. continuous function A(t, s) for s < t such that on its domain 

|o(t,s,y)F(e,s,y,z,e) - »(t,s,y*)F(0*, s,y*,z*, e)| 

< V 0 (€)4(t, s) (I e-0*| + I y-y*| + I z-z*| ) . 

There is a continuous function 5(t, s), s < t, such that on its 
domain 


| A(e(t-s),z)G(0,s,y,z,e) - A(e(t-s) , z*)G(0*, s,y*, z*, e) | 
< S(t, s) (| e-e*| + |y-y*| +|z-z*|). 


(5) There are positive constants M, Kp Kg and 

e, (0 < € - <e ) such that 

1 v 1 - o' 

t t 

/ A(t, s)exp / L(u, e^)du ds < Kp 
-oo s 

t t 

/ S(t, s)exp / L(u, e^du ds < Kg, t in R, 

-oo s 

•where L(u, e) = (2M+l)^(u, e) + (M+l)N. 

For' 0 < r ± < a*, 0 < r 2 < Og, T = (TpT 2 ) we define S r 
to be all v = (vpVg) where v^(0,t) is continuous from R^"*"^" 
into R m with | D(t) V-^0, t) | < fp where Vg(0,t) is continuous 
R^ +1 into R n with | E(t) Vg(0, t) j < Tg and where VpVg 


from 
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have period co in 0 and satisfy Lipschitz conditions ip 0 with 
constant M. For v in S we define H (0,t, e) 

= H(0,t,v^(0,t), v 2 (0,t),e) and F V , G V similarity. We assume for 
some x the following conditions hold: 

(6) There are continuous functions V-]_(£) on [0,e^], 

V^O) = V 2 (0) =0 so that . 

t 

/ |D(t)$(t, s^v-^Cs), s)F V (0(s),s,e)| ds < V^e), 

-00 

t 

e / |E(t)A(e(t-s), v 2 (0 (s) , s) ) G V (0 (s) , s, e) | ds < V 2 (e) 

-00 

for all v in S*" and continuous functions 0 from R into R m . 

(7) For any continuous function 0 from R into R m and 
v in S r let 

t 

P 0v (t) = / $(t, s,v 1 (0(s),s))F V (0(s),s,e)ds 
-00 

t 

qg v (t) = I A(e(t-s), v 2 (0(s),s))G V (0(s),s,e)ds. 

-00 

For any e in [0, e.^], fixed we assume 


|p ev « - W^i. iv<*> - 

are o(| t-t*| ) uniformly for v in S ' and 0 continuous from 
R into R m . 
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Theorem 5 . Assume that conditions (l) - (7) are satisfied. For e 
sufficiently small Ihere are functions (v 1 (6,t) ; v 2 (0,t) in S r 
such that y = v^(0,t), z = v 2 (0,t) is an integral manifold of ( 15 ). 

Proof : For v in let | V (t, t ; 0 o ) be the solution of 

(1*0 9 = h(e,t,v 2 (e,t)) + H v (e,t,€), 0 (t) = 0 O . 


The right side of the equation is Lipschitz in 0 with constant 
L(t, e) and is bounded by + C 2> A routine use of Gronwall's 
inequality gives 


(15) U v (t,T,e o ) - s T (t,t*,e*)| < 

\ 

T 

(|e 0 “ e 3 + ( C i + c 2 )|t-t*1)«p / L(u, e L ) du 

t 

for t > t and 0 < e < e^. Let k be a positive integer and let 
"be the functions (v^(e^t), v 2 (0,t)) mapping R - ^ x [-k,k] 
continuously into (R m ,R n ) with period oo in 0. Then S 1 " 
restricted to [-k,k] is a closed convex subset of B^. For v 
in S r define Tv(0,t) = (w^(0,t), w 2 (0,t)) by 

t 

w,(0,t) = / $(t, s,v (4 V (s,t,0),s)F V (£ V (s,t,0),s, e)ds 

-k 

t 

w ? (0,t) = e / A(e(t-s), v ? (| V (s,t,0), s)G V (| V (s,t,0), s,e)ds 
-k 
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for 0 in R^, -k < t < k. Tv exists by (6) since D,E are 
nonsingular. Since the right side of (l4) has period in 0 
we note £ V (s, t, 0 -kd) = £ V (s, t, 0) + to; consequently Tv has period 
a> in 0. For e sufficiently small v^e) < y^, v 2 (e) < r 2 
and hence | D(t) w^(0, t)| < y^, | E(t) w 2 (0, t) | < y 2 . We have 


Iw^e^t) - w^e^t*)! 

' t 

<V (e) / A(t,s)[2M+l]U v (s,t,e) - e (s,t*,e*)|ds 
° -k 


and a similar equation for w 2 (0,t). By (15) and condition (5) we 
have the family Ts"^ is equi continuous in (0, t) and for t = t* 
above for 6 sufficiently small |w^(0,t) - w^ (0,t)| < Mj 0-0*|, 
i = 1, 2. Note the condition on e is independent of k. 

Schauder's fixed point implies there is a fixed point of T, Let 
v be such a point and T > e 0 136 fixed. A repetition of the 
argument in Theorem 1 shows 0(t) = | v (t,T,0 Q ), y(t) = v (4 v (t, T,0 Q ),t), 
z(t) = v p (| v (t, t,0 ),t) is a solution of (13). 


For each positive k let v be a fixed point of T and 
let 0 k (t, x,0 o ), y^(t, t,©^), z k (t, t, 0 o ) "be the corresponding 
solutions of (13) . If Z is the set of positive integers, 
(v*(e,t): k in Z) is uniformly bounded and equicontinuous (by 
condition (7) on R - ^ X [-1,1]) so there is a Z^ C Z so that 

k in Z^} converges uniformly on R^ X [ - 1, 1] . In this 
way we obtain a decreasing sequence Z^ = [k^, k^ 2 , . . .} such that 
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(v k (0,t): k -in Z .} converges uniformly on R - ^ X Let 

k i = k ii» i = z * = ^ k i^i=r then k in z *} 

£ 

converges uniformly on sets of the form R x C, C compact in R. 

Let v(0,t) be the limit function. The corresponding sequence 
k k k 

0 ,y , z also converges uniformly on compact sets in R; thus the 
limit functions are solutions of (13). Consequently 

y = v-^Rjt), z = v 2 (0,t) 

is an integral manifold of ( 13 ) . 

The usual form for the 0 equation in (13) is 0 = t] + 
H(0, t, y, z,e) where rj is a constant. We include an example 
which may occur naturally (say by using Newton's equations) in a 
mathematical model and show how these equations can be reduced to 
the form ( 13 ) featuring an h term in the 0 equation which is 
not constant. Consider a weakly coupled nonlinear system 

x + £»? = -ex + e 2 x(t,x, x,y) 

y = a(t)y + b(t,y) + eF 1 (t,x,k,y) 

where the y equation for e = 0 was described in section 3, P 

4 

is a positive number and F^,X are continuous on R . For e = 0 
the first equation can be solved in terms of elliptic functions 
sn, dn and cn with modulus 1/ */2. Let x = p cn 0, 
x = -p 2 s/p sn0 dn 0, then equation above takes the form 
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e = vTp p - eH^t^p), 

V 

y = a(t)y + b(t,y) + eF^Q, t, y, p) , 
2 2 2 

p = -ep sn 0 dn 0 - 6 G-^G^t, y, p) . 


Let 

1 4K 2 2 

T = -^ / sn s dn s ds, 

.0 

® 2 2 

y(9, z) = z / (sn s dn s - r)ds 
0 

where 4 k is the real period of the elliptic functions, then under' 
the coordinate transformation given in the method of averaging 
([ 5], Chapters 14-17), p = z + ev(0,z) .the differential equations 
become 


0 = >/p z + €H 2 (0, t,y, z, e), 

y = a(t)y + b(t,y) + & 2 (e, t,y,z,t), 

2 

z = -erz + e G(0, t,y, z, e) . 

Relatively mild hypotheses on X, give the required smoothness 
conditions on H 2 , F 2 , G, etc. 
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6. CONCLUSION 

The use of the D(t) matrices in forming the underlying 
function spaces a) allows the asymptotic character of the solutions 
of the unperturbed system to he applied to the perturbed systems 
as in Case 2 of the example in section 3 and b) allows special time 
dependence of the perturbation to be taken into account as in Case 3 
of the section 3 example. In comparison to the version of Theorem 1 
given by May [ 11] our method of proof results in an improvement of 

the allowable size of the perturbation; however, our proof does not 

, 

reveal any periodic or almost periodic character and does not give 
unique solutions 

No Lipschitz conditions were used in section 4, thus no 
study was made of the differences y(t,a, a.^ - y(t,a, a*). (Here 
y(t,a, a^) is a solution furnished by the conclusion of Theorem 2 
with y^(a,o:, a^) = a^. ) If such conditions are imposed then 
asymptotic estimates of these differences can be made. Harbertson 
and Struble [8], (Theorem 2) and Hale [4], (Lemma 2.3) give such 

i 

estimates for integral manifolds. Even if Lipschitz conditions 
are imposed, our method of proof gives slightly larger bounds for 
the size of the perturbation at the expense of uniqueness. 

We remark here that the hypotheses on the 
A(v(e) (t-s),y)G(0, s,y, z) term in Harbertson and Struble [8], 
page 271, can include v(e) as a multiplier on the left. This 
will allow a larger class of examples. 
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